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ABSTRACT 

We present a detailed X-ray study of the intracluster medium (ICM) of the nearby, cool-core galaxy cluster 
Abell 478 (z = 0.088), based on Chandra and XMM-Newton observations. Using a wavelet smoothing hardness 
analysis, we derive detailed temperature maps of A478, revealing a surprising amount of temperature structure 
for an apparently well relaxed cluster. We find the broad band Chandra spectral fits yield temperatures which 
are significantly hotter than those obtained with XMM-Newton, but the Fe ionization temperature shows good 
agreement. We show that the temperature discrepancy is slightly reduced when comparing spectra from re- 
gions selected to enclose nearly isothermal gas. However, by simulating multi-temperature spectra and fitting 
them with a single temperature model, we find no significant difference between Chandra and XMM-Newton, 
indicating that non-isothermality cannot fully explain the discrepancy. 

We have discovered four hot spots located between 30-50 kpc from the cluster center, where the gas tem- 
perature is roughly a factor of 2 higher than in the surrounding material. We estimate the combined excess 
thermal energy present in these hot spots to be (3 ± 1) x 10 59 erg. The location of and amount of excess energy 
present in the hot spots are suggestive of a common origin within the cluster core, which hosts an active galactic 
nucleus. This cluster also possesses a pair of X-ray cavities coincident with weak radio lobes, as reported in 
a previous analysis, with an associated energy of less than 10% of the thermal excess in the hot spots. The 
presence of these hot spots could indicate strong-shock heating of the intracluster medium from the central 
radio source - one of the first such detections in a cool core cluster. Using the high resolution of Chandra, we 
probe the mass distribution in the core and find it to be characterized by a logarithmic slope of -0.35 ± 0.22, 
which is significantly flatter than an NFW cusp of -1 and consistent with recent strong lensing results for a 
number of clusters. 

Subject headings: galaxies: clusters: individual (A478) - cosmology: observation - galaxies: clusters: general 
- X-rays: galaxies: clusters 



1. INTRODUCTION 

Clusters of galaxies are critical sites for investigating 
the interaction between galaxies and their environment 
- particularly their impact upon the gaseous intracluster 
medium (ICM). There is substantial evidence to demon- 
strate that the ICM has been subjected to non-gravitational 
heating and/or cooling, so as to break the simple, self- 
similar scal ing of cluster propert i es with mass; e.g. the 
L-T x (e.g. lEdge & Steward 119911 fArnaud & Evrardl 119991 
Fairlevetal. 2000), M„„ - Tv lMohr Mathieso n. & Evrarfl 
1999) and isophotal size-temperature (Mohr & Evrard 1997) 
relations. Furthermore, it is apparent that the ICM is system- 
atically under-d ense and more extended in less massive ha- 
los (e.g. Ponman. Cannon. & Navarro 1999; Sanderson et al. 
| 2003|IOsmond & Ponmanll2004HAfshordi. Lin. & Sanderson! 
2005), and ther e is a increasing excess of entropy in the ga s in 
cooler systems (Ponman. Sande rson. & Finoguenovi2 003). 

While these observations can be explained by the im- 
pact of feedback a ssociated with galaxy formation (e.g. 
IVoit & Ponmanll2~003l) . the precise details of the mechanisms 
which mediate this interaction remain unknown. Of critical 
importance is the role of radiative cooling, which could fuel 
star formation and gradually deplete the reservoir of gas in 
the cluster core. Despite the inherently unstable nature of 
cooling by thermal bremsstrahlung, gas in the dense, undis- 
turbed cores of galax y clusters appears not to be cooling at the 
expected level (e.g. Ta mura et al.ll2001t iPeterson et alJl2003l 
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iKaastra et al. | l2004T). leading to greatly reduced mass depo - 
sition rates (iMakishima et al.ll200ll iBohringer et al.l 120021) . 
A plausible explanation for this behaviour is that thermal 
conduction acts to transfer signific ant amounts of energy 
from the outer region s of the ICM ( Kim & Naravan r2003t 
Voigt & Fabian 2004). How ever, this process cannot oper- 
ate at all ( Khosroshahi, Jones, & Ponman 2004) or effectively 
enough (e.g. Wise, McNamara, & Murray 2004) to heat the 
core in all cases, especially at lower temperatures. 

Since many "coolin g flow" clus t ers are know n to harbor 
central radio sources (Burnslll990t lEilekl 12004 . it is likely 
that active galactic nuclei (AGN) can deposit energy di- 
rectly into the I CM, and thus offset at least some of the 
cooling. Recently. ICroston. Hardcastle. & Birkm shaw (120051) 
have demonstrated that heating from low power radio galax- 
ies could account for the steepening of the L-Tx relation 
in X-ray bright galaxy groups. Although clusters are much 
more massive, the amount of energy output by powerful radio 
sources is certainly sufficient to heat the gas by the required 
amount (Chura zov et al.ll2~002l) : what is lacking is a clear un- 
derstanding of exactly how this energy is coupled to the ICM. 

A number of clusters have been shown to possess cavi- 
ties, where expanding radio lobes have displaced the ambi- 
ent medium and left a depression in the projecte d X-ray sur- 
face brightness (e.g. see the compilation of Birzan ^t alJ2004t 
and references therein). However, the coolest gas is generally 
found beside these lobes (Fabi anet al.l20 00: McNama ra et alJ 
120011) and their essentially subsonic expansion is not obvi- 
ously heating the ICM by a significant amount. In the mean- 
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time, the role of AGN driven bubbles in modifying the ICM 
is opening up as a ver y promising avenue of research (see 
Gardini & Rickei 2004, for a recent review). 

Alternatively, energy may be transferred to the ICM 
by means of acoustic wave s generated by turbulence 
(Fuiita. Suzuki. & Wada 2004), or viscous dissipation of 
sound waves originating fr om a centr a l radi o source 
jRuszkowski et alJ l2004albl) . iFabian et all J2003I) have re- 
cently discovered ripples in the X-ray emission from the 
Perseus Cluster, which point to a continuous energy output 
from the central radio source that can balance cooling within 
the innermost 50 kpc of the core. 

In this paper we examine in detail the properties of the ICM 
in a nearby galaxy cluster, using X-ray data from two dif- 
ferent telescopes. This apparently well relaxed system has 
a large cool core and its central galaxy hosts an active radio 
source, affording an excellent opportunity to study its inter- 
action with the ICM. We assume the following cosmologi- 
cal parameters: Hq = 70 kms" 1 Mpc" 1 , fl m = 0.3 = 0.7. 
Correspondingly, at our adopted redshift for A478 of 0.088, 
1"= 1.65 kpc. Throughout our spectral analysis we have used 
XSPE C 11.3.0, incorporati ng the default solar abundance ta- 
ble of lAnders & Grevessd G?89). All errors are la, unless 
otherwise stated. 

2. DATA REDUCTION 
2.1. Chandra X-ray Data Reduction 

A 42.4 ks Chandra observation of Abell 478 was made on 
January 29, 2001 using the Advanced CCD Imaging Spec- 
trometer (ACIS) in FAINT mode. The data reduction and 
analysis were performed with Ciao version 3.0.2 and CALDB 
version 2.26. This release of the Chandra CALDB calibration 
files was the first to incorporate a full correction for the degra- 
dation in quantum efficiency of the ACIS detectors. There- 
fore, no further steps were taken to correct for its effects. 

Three separate light curves were constructed for the SI, S3 
and the remaining CCDs combined, using the recommended 
energy and time binning criteria 1 . A small flare was iden- 
tified in the light curve from the SI CCD and another was 
found in the combined I2/I3/S2/S4 light curve: both were ex- 
cluded, leaving 40.7 ks of useful data. Only events with ASCA 
grades 0,2,3,4,6 were used, and bad columns and hot pixels 
were excluded. In addition, those events associated with cos- 
mic ray afterglows were identified and removed. A correction 
was applied to the level 1 events file to allow for the ACIS 
time-dependent gain variation, using the tool "corr_tgain" 2 , 
and a new level 2 events file was generated by reprocessing 
this modified level 1 events dataset. 

Since A478 fills much of the Chandra field of view, we use 
the Markevitch blank-sky datasets to estimate the background 
contribution 3 . To allow for small variations in the particle 
background level between the blank sky fields and the A478 
observation, we rescaled the effective exposure of the back- 
ground datasets according to the ratio of count rates in the 
particle-dominated 7-12 keV energy band for the SI CCD, 
which is furthest away from the cluster. To avoid the bias 
caused by the presence of contaminating point sources in the 
A478 dataset, we identified and excluded such features using 
the following iterative scheme. 

Images were extracted in the 0.5-2.0 keV band for the SI 

1 http://cxc.harvard.edu/contrib/maxim/bg. 

2 http://asc.harvard.edu/cont-soft/software/corr_tgain. 1 .0.html 

3 http://cxc.harvard.edu/contrib/maxim/acisbg. 



CCD in both the main (A478) and blank sky datasets. A back- 
ground image was created by smoothing the blank sky field 
with a Gaussian of width 1 arcmin, to filter out Poisson noise. 
The main image was then searched for sources with the ciao 
task WAVDETECT, using this background. The source re- 
gions found were then masked out of both the main and blank 
sky datasets, and the remaining counts in the 7-12 keV band 
were summed. A rescale factor was then determined as the 
ratio of the net count rate in the background dataset to that 
in the main dataset. The effective exposure time of the blank 
sky dataset was then multiplied by this rescale factor and the 
process repeated, until no new source regions were found. A 
total of 6 contaminating source regions were identified and 
excluded from the S 1 CCD data in this way. We found that 
the blank sky data were 8% higher than in the A478 obser- 
vation, consistent with the value of 10% quoted bv lSun et alJ 
d2003l hereafter lS03l) for the same observation. 

2.2. XMM X-ray Data Reduction 

XMM-Newton observed A478 in orbit 401 (observation ID 
0109880101), using the thin filter to block visible light. In 
the analysis that follows we present the data obtained with the 
pn-CCD camera dStruder & et alJ20 01 ). operated in extended 
full-frame mode, with an exposure time of 42 ks. The ini- 
tial stages of the data reduction were performed using XMM- 
SAS 5.4.1. Screening of the data for contamination from 
flares was achieved by examining a light curve of detector 
counts in the 10-15 keV band. We have adopted the approach 
ba sed on the analysis o f the count rate histogram, as described 
in iZhang et alJ i2004f) . which is more sensitive to the back- 
ground conditions during the observations than using a fixed 
threshold to identify good time intervals. 

With these screened photon event files we produced a pn 
image of A478 in the energy bands 0.5-2.0 and 2.0-7.5 keV. 
In the imaging analysis, we included photons near the pn- 
CCD borders, near bad pixels and offset columns, in order 
to reduce the width of the gaps between CCD chips. This is 
a reasonable approach to adopt when dealing with qualitative 
images in broad energy bands. For the spectral fitting, de- 
scribed in 5j3j we exclude all these border events, since a small 
fraction of those photons have an incorrect energy, typically 
due to the registration of a double event pattern as a single 
event. After subtracting the expected out-of-time events from 
the image, we corrected the resid ual for vignetting and ex- 
posure using the latest calibrations (Lumb et al. 2003), which 
have been incorporated within the XMMSAS 6.0 release. 

Emission from A478 fills the entire XMM-Newton field of 
view and so it was not possible to determine the background 
level reliably from this dataset. Therefore, an observation of 
the Chandra Deep Field South was used to provide an esti- 
mate of the background. We have confirmed that this choice 
of dataset is well matched to the detector background in the 
10-15 keV band for our observation of A478. However, for 
the goal of the current analysis the choice of the background is 
less important, since the measured flux is strongly dominated 
by emission from the cluster. 

3. CHANDRA SPECTRAL ANALYSIS 

A global spectrum of A478 was extracted from the emission 
filling the entire S3 chip, excluding contaminating sources. 
The spectrum was grouped to a minimum of 20 counts per 
bin, to enable us to use the x 2 fit statistic. The CCD detector 
responses vary as a function of position within our spectral 
extraction region. We therefore generated composite response 
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FIG. 1. — Projected gas temperature as a function of radius. The diamonds 
are the XMM-Newton pn data and the barred crosses are the Chandra data. 
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FIG. 2. — Projected gas metallicity as a function of radius. The diamonds 
are the XMM-Newton pn data and the barred crosses are the Chandra data. 



files using the CIAO tasks MKWARF and MKRMF, by aver- 
aging the contributions made from those parts of the detector 
with different responses, weighting each component accord- 
ing to the distribution of counts in the 0.5-2.0 keV band. 

Using XSPEC 11.3.0, a single temperature MEKAL hot 
plasma model combined with a WABS galactic absorption 
component, was fitted to the data between 0.7 and 7.0 keV, 
to maximize the signal-to-noise ratio (S/N). This yielded a 
best-fit temperature of 6.56 ± 0.06 keV, with a metallicity of 
0.33 ±0.01 solar, and an absorbing column of (2.91 ±0.01) x 
10 21 cm" 2 . The redshift was also left free to vary, giving a 
best-fit value of 0.0879^j] w hich agrees well with the op- 
tical value of 0.0881 ± 0.0009 (IStruble & Roodlll999l) . We 
note that S03 found the following spectral parameters based 
on their 0.7-8.0 keV analysis of data from the S2, S3 & S4 
CCDs combined using this observation: kT = 7.18 ±0.11 
and jVh = (2.59 ±0.03) x 10 21 cm -2 , with an abundance of 
0.37 ± 0.02 solar (all at 90% confidence). The differences 
are not surprising, since the gas temperature shows a marked 
decrease in the cluster center and the metallicity and absorb- 
ing column both exhibit central enhancements (see £13- li be- 
yond the S3 chip, to which we restrict our analysis. However, 
the significant discrepancy in absorbing column is mainly 
due to the improvement in the Chandra calibration since the 
ISun et all analysis - a factor which we will return to in Sec- 
tional 

3.1. Chandra vs XMM 

In view of our intention to exploit the complementary prop- 
erties of Chandra and XMM-Newton, by combining spectral 
data from both (see |4j, we present here a detailed comparison 
of results obtained from the two observatories. We are further 
motivated by the apparent discrepancy between the tempera- 
tures inferred from Chan dra and XMM-Newton obse rvations 
of A478, as reported by Pointecouteau et al. (2004). Given 
the potential for spatial variation in cluster gas properties (as 
well as in the absorbing column), we perform a direct com- 
parison of spectral properties in a series of contiguous annuli. 



Ten logarithmically-spaced radial bins were defined out to a 
maximum radius of 12' (see Tabled, centered on the peak of 
the emission as measured using Chandra (RA. 04 h 13 m 25 s .2, 
decl. 10°27'53"). The innermost bin radius was increased 
to 21", to reduce the effect of scattering between bins, due to 
the comparatively broad point spread function (PSF) of XMM- 
Newton. Equivalent Chandra spectra were extracted in the 7 
innermost annuli fully covered by the S3 chip. 

Absorbed MEKAL models were fitted to each spectrum, 
as before, and radial profiles of temperature, metallicity and 
absorbing column are plotted in Figs.^J3J respectively, show- 
ing the Chandra data as barred crosses and the XMM-Newton 
data as diamonds. The temperature profile is typical of a 
'cool-core' cluster, showing a sharp drop towards the center, 
within The region outside the core is roughly isothermal, 
with some indication of a negative gradient further out. It is 
also clear that the Chandra points are systematically higher 
in temperature compared to XMM-Newton. This discrepancy 
persists in the measured abundances, with Chandra again giv- 
ing somewhat higher values than XMM-Newton in all annuli 
(Fig.|2ji. Nonetheless, both observations show an increase in 
the gas metallicity in the cluster core, approximately coinci- 
dent with the observe d decrease in temperature in this region. 
Vikhli nin et alJ J2004I) present a recent analysis of the A478 
Chandra data, which incorporates the latest calibration. Their 
projected temperature profile is similar to our own, although 
it extends beyond the S3 chip, reaching a slightly hotter peak 
of ~9 keV at roughly 150 arcsec. 

Furthermore, the inferred absorption column also rises to- 
wards the center, and there is close agreement between the 
Chandra and XMM-Newton results (Fig. |3J. Also plotted is 
the galactic absorpti on of 1.51 x 10 2 1 cm" 2 based on HI ob- 
servations (Dickev & Lockman 1990) interpolated to the cen- 
ter of A478. It can be seen that the fitted points lie well 
above this line, implying an observed absorption roughly 
twice as high. However, this excess is not localized to the 
cluster core, and is probably associ ated with absorption in our 
galaxy JPointecouteau et afll2004). This behavior is not un- 
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TABLE 1 

Spectral fitting results (0.7-7.0 keV) 



Rout" (") 




Chandra 






XMM-Newton 






kT (keV) 


Z (solar) 


nH{l0 21 cm 4 ) 


X 2 /dof 


kT (keV) 


Z (solar) 


nH (10 21 cm 4 ) 


X 2 /dof 


21 


4.25 ±0.08 


0.47 ±0.03 


2.98 ±0.06 


459.8 / 346 


4.43 ±0.06 


0.32±0.01 


2.87 ±0.04 


1367.8/ 1006 


43 


6.14±0.16 


0.35 ±0.03 


2.86 ±0.05 


434.7 / 376 


5.67 ±0.07 


0.28 ±0.01 


2.83 ±0.04 


1456.2/ 1150 


61 


6.76 ±0.20 


0.36 ±0.04 


2.81 ±0.06 


415.3/357 


6.19±0.12 


0.24 ±0.01 


2.76 ±0.04 


1312.2/1083 


87 


7.47 ±0.20 


0.29 ±0.03 


2.76 ±0.06 


437.6 / 374 


6.66 ±0.13 


0.23 ±0.01 


2.71 ±0.04 


1197.7/ 1095 


123 


7.67 ±0.28 


0.30 ±0.03 


2.74 ±0.06 


455.3/379 


6.57 ±0.12 


0.26 ±0.02 


2.69 ±0.04 


1325.2/ 1096 


175 


7.91 ±0.31 


0.26 ±0.04 


2.64 ±0.06 


443.5 / 385 


7.06±0.13 


0.23 ±0.02 


2.62 ±0.04 


1189.4/1116 


250 


8.12±0.37 


0.27 ±0.05 


2.67 ±0.07 


454.1 /379 


6.93 ±0.16 


0.20 ±0.02 


2.57 ±0.05 


1039.3/ 1058 


355 










6.96 ±0.21 


0.25 ±0.03 


2.20 ±0.06 


1025.5 / 1002 


506 










6.87 ±0.35 


0.32 ±0.05 


2.18 ±0.11 


1095.8/ 975 


720 










5.92 ±0.57 


0.38±0.11 


2.40 ±0.25 


1035.0/ 1030 



NOTE. — "Outer radius of annulus. The redshift was fixed at 0.088 for all spectra. All errors are 1<t. 
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Radius (arcsec) FIG. 4.— The 6.0-6.8 keV ACIS spectrum, best fitting MEKAL model 

and residuals from annuli 4-7 combined (see text for details). 

FIG. 3. — Projected absorbing column as a function of radius. The dia- 
monds are the XMM-Newton pn data and the barred crosses are the Chandra 
data. The HI inferred absorption value is indicated by the dotted line. 



expected, since HI measurements become an increasingly un- 
reliable predictor of the total galactic absorption for columns 
exceeding ~6x 10 20 cm" 2 due to contributions from molecu- 
lar hydrogen, for example (Lockman 2004). An excess ab- 
sorpti on was also i nferred from photometric observations by 
Garilli et al. ( 1996), who noted that the early-type galaxies in 
A478 are much redder than expected. 

In light of the disagreement in both temperature and metal- 
licity between Chandra and XMM-Newton, we performed ad- 
ditional tests to verify the impact of any calibration errors on 
our results. Specifically, we checked the ionization temper- 
ature of the ICM, as measured from the prominent iron line 
redshifted to ~6.1 keV in the observer frame. We extracted a 
spectrum from annuli 4-7 combined, thus excluding the inner 
cluster core, which contains gas at a wide range of tempera- 
tures and metallicities, as already seen. This spectrum was fit- 
ted with an absorbed MEKAL model as before, in the ranges 
0.7-7.0 keV and 6.0-6.8 keV (see Fig.E). Since the absorb- 



ing column was poorly constrained when fitting the narrow 
(high) energy range, we fixed it at the value obtained from the 
0.7-7.0 keV fit. However, we also tried fixing Nh at the galac- 
tic value deduced from HI observations (1.51 x 10 21 cm" 2 ). 
This made no appreciable difference to the results, which are 
summarized in the left half of Table |2] The corresponding \~ 
and degrees of freedom are (71 1.07/427) and (39.49/51), for 
the broad and narrow range Chandra fits, respectively. It is 
clear from the broad band fit that, even within this apparently 
roughly isothermal region, the cluster emission is not well de- 
scribed by a single temperature plasma. This hints at the pres- 
ence of significant spatial variation in the temperat ure o f the 
ICM - a possibility which we investigate further in S15.3I 

For direct comparison, the equivalent test was performed on 
the XMM-Newton pn spectrum from annuli 4-7; the results 
are summarized in the right hand side of Table [2] As with 
Chandra, it was necessary to fix Nh for the narrow band fit at 
the optimum value obtained from the 0.7-7.0 keV fit. Once 
again we checked that the results do not change noticeably 
when Nh was fixed at the HI value. However, we did find that 
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TABLE 2 

Spectral fitting results for annuli 4-7 combined 







Chandra 








XMM-Newton 






Fit Range (keV) 


kT (keV) 


Z (solar) 


rcff (10 21 crrr 2 ) 


X 2 /dof 


kT (keV) 


Z (solar) 


nH{W 21 cm" 2 ) 


X 2 /dof 


0.7-7.0 


7.73±0.15 


0.28 ±0.02 


2.72 ±0.03 


711.1 /427 


6.81 ±0.06 


0.28 ±0.01 


2.66 ±0.02 


1338.0/ 1256 


6.0-6.8 


6.49 ±0.61 


0.22 ±0.05 


2.72 (frozen) 


39.5/51 


6.69 ±0.26 


0.24 ±0.02 


2.66 (frozen) 


158.5/152 



NOTE. — The redshift was fixed at 0.088 for Chandra, but left free to vary for the XMM-Newton, broad-band fit; a simple linear gain offset was 
fitted to the narrow-band XMM-Newton spectrum to give the correct redshift (see text for details). All errors are Icr. 



freezing the redshift at the optical value of 0.088 produced an 
unacceptable fit (x 2 /dof = 423.2 0/154), indicative of a small 
cali bration offset , also found by iPointecouteau et a"fl (120041) 
and de Plaa et al1d2004l) . This is due to the uncertainty in the 
current CTI correction affecting the pn detector in extended 
full frame mode only (K. Dennerl 2004, private communica- 
tion), a correction for which has only just been incorporated 
into the XMMSAS 6.1 release. Leaving the redshift free to 
vary produced a best-fit value of 0.074 ±0.001. To provide 
a fairer comparison with Chandra in a region heavily domi- 
nated by emission line flux, we therefore applied a small mod- 
ification to the response matrix file to improve the fit, using 
the XSPEC "gain" command. This produced an acceptable fit 
with a x 2 /dof of (158.65/154), compared to (1337.99/1261) 
for the broad range fit, without the gain modification. The 2 
fitted gain parameters were frozen at their best fit values in the 
calculation of MEKAL parameter errors. 

Table |2 shows that the broad-band Chandra spectrum pro- 
duces an anomalously high temperature, which disagrees with 
the Chandra ionization temperature measurement, as well as 
with both the broad-band and ionization temperature fits from 
XMM-Newton. This could indicate a problem with the ACIS 
calibration as the source of the discrepancy between Chan- 
dra and XMM-Newton reported by Pointecouteau et al. ( 2004) 
and seen in Table [2] for this dataset. However, this behavior 
can be explained by significant temperature variation within 
each annulus, giving rise to different characteristic average 
temperatures when folded through the Chandra and XMM- 
Newt on re sponses and fit with a single phase plasma model 
(see 35. 3> . The implication is that better agreement between 
Chandra and XMM-Newton ought to occur for spectra ex- 
tracted from appro xima tely isothermal regions, which we find 
to be the case (see £15 .31 . 

In the meantime, we have neglected to correct for this bias 
in what follows, in view of the impracticality of restricting 
our Chandra analysis to a small, high energy range, and given 
the somewhat better agreement between Chandra and XMM- 
Newton seen in the deprojected temperature profile (see Sj4}. 

3.2. The Effects of Multiphase Gas 

In order to shed some light on the discrepancy between 
Chandra and XMM-Newton described above, we have con- 
ducted a series of simulations of multiphase spectra, which we 
have fitted with a single temperature spectrum. We have used 
the spectral responses and background spectrum from annuli 
4-7 (see Tabled) for both telescopes and assumed a value 
of galactic absorption of 2.7 x 10 21 cm" 2 for all these simu- 
lations. Initially we generated 1000 realisations of a spec- 
trum comprising 4 temperature components (6,6.5,7.5 & 8 
keV, weighted in the ratio 7:10:30:7), with an abundance of 
0.3 solar. The results of fitting this spectrum with a single ab- 



sorbed MEKAL model recovered the following mean temper- 
atures and standard deviations: 7.1 1 ±0.10 (Chandra 0.7-7.0 
keV); 7. 12 ±0.06 (XMM 0.7-7.0 keV); 7. 16 ±0.45 (Chandra 
6.0-6.8 keV); 7. 17 ±0.19 (XMM 6.0-6.8 keV). The excellent 
agreement between these 4 cases demonstrates that moder- 
ately multiphase gas with the same metallicity cannot account 
for the discrepancies seen in Table |2] despite the differences 
in the spectral responses for Chandra and XMM-Newton. 

We then simulated a spectrum with 2 temperature compo- 
nents (5 and 12 keV) and an abundance of 0.3 solar (case A), 
0.2 & 0.5 solar (case B), and 0.5 & 0.2 solar (case C), respec- 
tively. The fitting results are summarised in Table|3] It can be 
seen that there is consistently good agreement between Chan- 
dra and XMM-Newton in all 3 cases. Furthermore, it can also 
be seen that the 0.7-7.0 and 6.0-6.8 results are generally in 
good agreement when both have the same metallicity. How- 
ever, when the hotter phase has higher metallicity, the narrow- 
band ionization balance fit yields a significantly hotter tem- 
perature for both Chandra and XMM-Newton. Conversely, 
when the cooler phase is more enriched, the 6.0-6.8 keV fit 
yields a cooler temperature, albeit at rather lower significance. 

As an aside we note that, in all cases, the metallicity ob- 
tained with a single-temperature fit to the multiphase spectra 
overestimates the input values, by ^25% in the simple case 
where the 2 phase have the same abundance. Where the metal- 
licities of the 2 phases differ, the abundance is over-estimated 
by an even greater amount, for both the broad- and narrow- 
band fits. 

Based on these simulations, we therefore conclude that the 
discrepancy we observe between Chandra and XMM-Newton 
cannot readily be explained by the presence of multiple tem- 
perature components, even with differing abundances. 

4. DEPROJECTION ANALYSIS 

To determine the intracluster gas properties, we imple- 
mented a standard "onion peeling" scheme to deproject the 
X-ray data, under the assumption of spherical symmetry. 
The spectral deprojection was performed using the PROJCT 
model in XSPEC, using a single-temperature MEKAL plasma 
model with a galactic absorption component (WABS). Only 
data within 0.7 and 7.0 keV were used, and each spectrum 
was grouped to a minimum of 20 counts per bin, including 
the background contribution. The foreground absorbing col- 
umn density was left free to vary in each annulus, as was the 
MEKAL temperature, abundance and normalization. 

The XSPEC MEKAL normalization is defined as 

10" 14 [ 

K= a 7n 7uL ^2 WndV, (1) 
4ir(D A (\+z)Y J 

where Da is the angular diameter distance, z is the redshift 
and dV is the volume from which the deprojected emission 
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TABLE 3 

Spectral fitting results for simulated multiphase spectra 



Sim. 


Fit Range (keV) 




Chandra 








XMM-Newton 




kT (keV) 


Z (solar) 


nH(W 21 cnr 2 ) 


X 2 /dof 


kT (keV) 


Z (solar) 


nH (10 21 cm" 2 ) 


X 2 /dof 


A 


0.7-7.0 


7.28±0.12 


0.37 ±0.02 


2.66 ±0.03 


439.1 /427 


7.36 ±0.06 


0.38 ±0.01 


2.65 ±0.02 


1294.6/ 1260 


A 


6.0-6.8 


7.65 ±0.63 


0.39 ±0.08 


2.70 (frozen) 


51.8/51 


7.70 ±0.28 


0.39 ±0.04 


2.70 (frozen) 


157.7 / 155 


B 


0.7-7.0 


7.50±0.15 


0.38 ±0.03 


2.65 ±0.04 


443.4/427 


7.67 ±0.09 


0.39 ±0.01 


2.61 ±0.02 


1324.4/1260 


B 


6.0-6.8 


9.19 ±0.57 


0.54±0.10 


2.70 (frozen) 


51.5/51 


9.16±0.22 


0.53 ±0.04 


2.70 (frozen) 


156.9 / 155 


C 


0.7-7.0 


7.03 ±0.11 


0.49 ± 0.03 


2.67 ± 0.03 


455.7 / 427 


7.07 ±0.06 


0.48 ±0.01 


2.67 ±0.02 


1331.4/ 1260 


C 


6.0-6.8 


6.45 ±0.53 


0.41 ±0.08 


2.70 (frozen) 


51.5/51 


6.46 ±0.24 


0.39 ±0.04 


2.70 (frozen) 


157.6/155 



NOTE. — Input spectra comprise equally-weighted 5 & 12 keV components with equal abundance of 0.3 solar (A), 0.2 & 0.5 solar (B) and vice versa (C), 
respectively. Errors are the standard deviation of the 1000 Monte Carlo realisations in each case. 
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FIG. 5. — Gas electron number density as a function of radius. Errors on 
the density are too small to be plotted and the radial bin width error bars have 
been omitted for clarity. The hollow circles represent the Chandra points 
and the XMM-Newt on data are the "+ " symbols. Also shown is the ROSAT 
density profile from Mohr et al. 1 1999), plotted as a solid line for clarity. 
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FIG. 6. — Deprojected temperature as a function of radius. The diamonds 
represent the XMM-Newton pn data and the barred crosses are the Chandra 
data. The dotted crosses are from the S03 analysis of the same Chandra ob- 
servation, using older calibration data. Although some discrepancy between 
Chandra and XMM-Newton remains, it has been reduced with improvements 
to the CALDB and is less pronounced than in the projected T(r) fFig.lTi. 



originates. The best-fit normalization for each spectrum was 
converted directly into a mean gas density in the correspond- 
ing spherical shell, assuming a number density ratio given by 



«e = «H = 1.20 «H, 

Me 



(2) 



appropriate for a fully ionized plasma with a metallicity of 0.3 
times the solar value, using our adopted abundance table from 
lAnders & Grevessel(ll989i) . 

Data from the outermost bin were discarded, since this vol- 
ume incorporates a contribution from emission beyond the 
last shell that implies a non-trivial geometry for the pur- 
poses of determining the mean gas density. This leaves a 
total of 6 Chandra and 9 XMM-Newton radial bins, span- 
ning the innermost 8.4 ' (835 kpc) of the ICM. Following 
iLewis. Buote. & Stockd d2003l) we assign an effective radius 
to each annulus, which is approximately equivalen t to an 
emission-weighted mean, given by (see McLaughlin 1999) 



0.5 



I ' 01 



/2 + r 3/2 
out ^'in 



2/3 



(3) 



The resulting density profile is plotted in Fig. |5] a nd also 
includ es the data from the ROSAT analysis of Mo hr et ail 
( 1999), shown as the solid line. There is generally good 
agreement between the 3 sets of data, although the Chandra 
points appear to be systematically higher than XMM-Newton. 
It is particularly encouraging to note the close agreement be- 
tween the ChandralXMM-Newton values and the ROSAT pro- 
file, since the latter was derived from an analytical (double 
/3-model with PSF corrections) fit, as compared to the non- 
parametric "onion peeling" method used in this analysis. 

The corresponding deprojected gas temperature profile is 
shown in Fig. |6] with the Chandra and XMM-Newton points 
plotted as before. It can be seen that the Chandra tempera- 
tures are hotter, except for the innermost bin, where PSF scat- 
tering acts to smooth out the gradient in the XMM-Newton 
data. However, the agreement between the 2 observations is 
better than in the projected kT(r), with an overlap in the lcr 
error bounds for all but 3 annuli. The deprojected temperature 
gradient is steeper than that observed in the projected profile 
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(Fig. [Q, as expected (since projection smooths out any such 
gradients), and reaches a minimum central valu e of ~ 3.5 ke V. 
Also shown, for comparison, are the points from S03 using the 
same Chandra data, plotted as dotted crosses. These temper- 
ature values are systematically hotter than our measurements. 
We attribute this discrepancy to changes in the Chandra cali- 
bration between CALDB 2.15 and CALDB 2.26, particularly 
with respect to the method for handling the effects of the 
ACIS contamination, which causes an energy-dependent ab- 
sorption that can affect temperature estimates. As can be seen, 
the newer ACIS calibration produces temperatures which are 
more consistent with those from XMM-Newton. 

5. SPECTRAL MAPPING ANALYSIS 
5.1. Hardness Method 

Although the surface brightness of A478 appears to exhibit 
a smooth morphology typical of a relaxed cluster, it is difficult 
to gauge the true state of the gas from studying images alone. 
Here we describe the process of creating a temperature map 
from the Chandra X- ray data, based on a wav elet-smoothed 
hardness ratio image (Finoguenov et al. 2004b). The key ad- 
vantage of this approach is that it allows fine scale tempera- 
ture structure to be resolved, without any prior knowledge of 
its morphology, since we are able to work with high resolution 
images of the cluster. 

We select two hardness bands, with the following energy 
ranges: 0.5-2.0 keV (soft band) and 2.0-7.0 keV (hard band). 
The hardness ratio (given by the hard band flux divided by 
the soft band flux) depends strongly on the gas temperature 
and the absorbing column along the line of sight, as well as 
more weakly on the metallicity of the gas. Under the assump- 
tion that the variation in absorbing column over the S3 chip 
is small (see the barred points in Fig. [3}, we can therefore use 
the hardness ratio to estimate the temperature directly. 

Taking the best-fit model to the global spectrum (§3), we 
fix the abundance and absorbing column and generate a series 
of 200 new models with different temperatures in the range 
0.5-15.0 keV. For each of these spectra, we evaluate the ra- 
tio of the predicted model flux in the hard to soft bands. We 
then fit a 6th order polynomial curve to the data, to obtain 
an expression for the gas temperature as a function of hard- 
ness. The residuals from this best-fit relation are less than 1% 
throughout the range of hardness values to which we apply it. 

To create a hardness ratio map from X-ray data it is nec- 
essary to apply some smoothing to suppress Poisson fluctu- 
ations. This ensures non-zero pixel values in the denomi- 
nator (soft band image) but can also reveal subtle variations 
in hardness on small scales. We use a wavelet decompo- 
sition approach, which is ideal for highlighting low surface 
brightness diffuse emission, even in the presence of bright 
embedded point sources, whilst retain ing complex structural 
information (see Vikhlin in et al.lll99 8. for details). This ap- 
proach also provides information on the significance of struc- 
tures identified on each scale and has recently been used in 
the analysis of XMM-Newton observations of several clus- 
ters (Finoguenov et al. 2004b a; Henrv, Finoguenov, & Briel 
120041 

Initially we extract source and background images and cor- 
responding exposure maps in both the soft and hard bands for 
the whole S3 chip. The corresponding background image is 
then subtracted from each hardness image, and the resulting 
image divided by the appropriate exposure map. The expo- 
sure maps are obtained by weighting the contributions to the 
effective area of the detector by the best-fitting model to the 



global spectrum fitted in |J3] within each of the soft and hard 
bands separately. Pixels with an equivalent exposure less than 
25% of the highest value within the S3 chip were masked out 
to improve the S/N of the resulting hardness image. 

A wavelet transform decomposition was then applied sep- 
arately to the soft and hard band images, over 6 different 
scales sizes, increasing in integer powers of 2 from 2 2 pix- 
els to 2 7 pixels, corresponding to angular sizes of ^2" and 
respectively. A 5a threshold was used to identify gen- 
uine sources, whose extent was defined to enclose those 
neighborin g pixels lying above a 2a threshold. Following 
Finog uenov et alJ (|T004b), an additional smoothing was ap- 
plied to each wavelet transform in turn, by convolving it with 
a Gaussian kernel of width equal to that scale, to reduce dis- 
continuity artifacts associated with separating the image into 
different scales. The soft and hard band images were recon- 
structed by summing all six smoothed wavelet transformed 
images obtained at each of the scale sizes. The final hardness 
map was obtained as the ratio of the hard to soft images. Each 
pixel was then converted into a temperature value by using the 
polynomial relation derived above, with unphysical hardness 
values (confined to a few pixels near the S3 chip boundary) 
masked out of the resulting temperature map (i.e. values ly- 
ing outside the fitting range of the calibration curve). 

5.2. Spectral Fitting Method 

There are several key assumptions of the hardness-based 
approach to generating a temperature map that could produce 
incorrect results. Specifically, there may be subtle, small- 
scale variations in absorbing column and gas metallicity that 
could manifest themselves as spurious features in the inferred 
temperature distribution. Here we present a more sophisti- 
cated method to determine spatial variations in gas tempera- 
ture, in order to assess the validity of the above technique. 

With an estimated temperature map we are well positioned 
to identify interesting features in the image for targeting with 
a conventional spectral fitting approach. Using contour maps 
of the temperature and soft band images, we defined a series 
of 53 separate regions, each containing a minimum of 1000 
counts in the raw image, originating from gas of similar pro- 
jected temperature and surface brightness, as implied by our 
wavelet decomposition analysis. 

We extract source and background spectra and response 
files in each of these 53 regions, and fit each with an ab- 
sorbed MEKAL model as before. Images of these regions, 
showing the best-fit temperature and la bounds, are plotted 
in the lower half of Fig.0 The color table for the spectral fit 
temperature map is identical to that used in the hardness tem- 
perature maps in the upper half of the figure. The broad fea- 
tures in the hardness map are confirmed in the spectral map: 
specifically, the temperatures of the cool core and surrounding 
regions agree very well. Moreover, the spectral analysis con- 
firms the existence of 4 prominent hot spots located around 
the cluster center, with temperatures as high as 12 keV - al- 
most double that of the surrounding medium in all cases, al- 
beit with large uncertainties. The origin of these hot spots is 
unclear, but it seems unlikely that they can all be attributed 
to point source contamination: visual inspection of the raw 
image reveals no evidence for a surface brightness excess at 
these locati ons. We will return to the interpretation of these 
features in 

5.3. Spectral Mapping Deprojection 
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FIG. 7. — Wavelet smoothed temperature maps of A478 from Chandra (top left) and XMM-Newton (top right), derived from the hardness ratio method. 
The positions of 4 hot spot regions are indicated by the letters A-D, and the cavities discovered by S03 are also labeled. The bottom half shows the Chandra 
temperature map derived from the spectral fitting method; from left to right the panels show the best-fit values, la low and upper limits on the temperature. The 
innermost green/yellow/orange regions are the hot spots. 



Using our spectral mapping results, it is possible to per- 
form an approximate deprojection, to provide a comparison 
with the "onion-peeling" method. To achieve this, we assume 
that all the emission projected onto each region originates in 
a spherical shell bounded by the inner- and outermost radii 
enclosing the region, measured from the surface brightness 
centroid used previously. The volume used to calculate the 
gas density from the MEKAL normalization (equation 0, is 
formed from the intersection of this shell with the cylinder 
having the cross section of that region. A more detailed de - 
scription of this approach is presented in Henry et al. (2004). 

In this "deprojection" scheme, no attempt is made to model 
out or otherwise remove the contaminating emission from 
regions outside of the volume element under consideration, 
which are projected onto the spectral region. Therefore, we 
take the projected temperature and metallicity information to 
apply to the 3 -dimensional volume element. However, in se- 
lecting these regions, we have ensured that that gas is essen- 
tially isothermal as measured in projection, which minimizes 
the additional smoothing effect of averaging over a range of 
temperatures. Correspondingly the inferred temperatures and 
abundances are more likely to be representative of the depro- 
jected gas properties. 

A plot of the temperature profile from the mapping analysis 
is presented in Fig. [8] with the radial axis expressed in units 
of R500 ( 1 349 kpc ), as determined from our mass model of the 
cluster (see ^15 .4-i . It can be seen that there is a large degree 
of temperature scatter, with the hot spots clearly visible, al- 
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FIG. 8. — Temperature profile from the spectral mapping analysis. No dis- 
crepancy between the Chandra (diamonds) and XMM-Newton (solid circles) 
points is evident in the mapping data. 



beit with large errors. The overlap between the Chandra and 
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XMM-Newton spectral mapping data is somewhat limited, due 
to the poorer spatial resolution of XMM-Newton combined 
with the impact of edge effects in the wavelet reconstruction 
limiting the outer radius of data from the Chandra S3 chip. 
However, it can be seen that the discrepancy between Chan- 
dra and XMM-Newton found in the annular spectral fits is not 
apparent in the spectral mapping data. Although this is due 
in part to the larger uncertaintied in the data points, it may 
also point to the source of the discrepancy being the varied 
temperature structure in the ICM, which gives rise to a range 
of temperatures within a given annulus. Since XMM-Newton 
is more sensitive to cooler temperatures and Chandra is more 
sensitive to hotter ones, such a situation is likely to lead to 
a bias in fitting a single temperature spectrum, as observed. 
Moreover, the iron ionization fit is more sensitive to cooler 
gas phases, which emit proportionately more line flux, which 
explains why the Chandra 6.0-6.8 keV result (Table |2ji was 
in good agreement with XMM-Newton. 

Fig.|9]shows the entropy profile from the spectral mapping 
analysis. A lso shown, as a dashed line, is a fit of the form 
S oc r l 1 (c.f. Tozzi, Scharf, & Norman 2000), which has been 
shown to provide a good mat ch to the outer regio ns of most 
reasonably relaxed clusters (iPonman et al.l 12003). To high- 
light the trend in the data, the points have been smoothed us- 
ing a locally weighted regression in log-log space (solid line). 
This technique smooths the data using a quadratic function 
which is moved along the set of points to build up a curve, in 
an analogous fashion to how a moving average is computed 
for a time series. The algorithm used is implemented in the 
LOWES S function in th e R Project statistical environmen t 
package 4 (version 2.0.0) (R Developme nt Core Tearr] l2004). 
and further details can be found there. The agreement in the 
outer regions between this curve and the S oc r l 1 line is strik- 
ing, although the data deviate slightly above this simple rela- 
tion in the cool core (within ~0. 1 R500 - see Fig.|8j. We defer 
further discussion of this issue to j 



5.4. Mass Distribution 

The deprojected gas T(r) and p(r) can be used to infer the 
gravitating mass profile, assuming hydrostatic equilibrium, 
given by 



M„ rav (r) : 



kT(r)r 
Gpm p 



d In p d In T 
+ 



dlnr dlnr 



(4) 



(e.g. iFabricant. Lecar. & G orenstein 1980), where p is the 
mean molecular weight of the gas and m p is the proton mass. 
In order to evaluate analytically the gradients in this equation, 
we fit a 3rd order polynomial to both the T(r) and p(r) data in 
log space. The resulting mass profile is evaluated at the radii 
of the input T(r) and p(r) values and is well described by a 
profile of the form 



p(r) = 



x n (l+x) 



3-n • 



(5) 



(e.g. lZhaollT996l) . where po is the central density, x = r/r s , 
and r s is a characteristic scale radius. For r -C r s the 
profile is characterized by a cusp, with p oc r~". Equa- 
tion |5]isJhe_gejierahzedform^ the so-called NFW profile 
( Nava rro. Frenk. & WhitdU995l) . which is obtained for n = 1 . 
A reasonable proxy for the virial radius of the halo can be ob- 
tained from r2oo, the radius enclosing a mean overdensity of 

4 http://www.r-project.org 
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FIG. 9. — Entropy profile from the spectral mapping deprojection. The 
point styles are the same as in Fig.|S] The dashed line is an unweighted fit to 
the data, of the form S oc R L1 . For comparison, the solid line shows a locally 
we ighte d regression of the data in log-log space, using the method described 



200 with respect to the critical density, p C rit(z)> of the Universe 
at the observed cluster redshift, z, given by 



Pcrit(z) = E(z) 



2 3// o 
8vrG' 



where G is the gravitational constant, Ho is the Hubble con- 
stant at z = and 



£( z ) = (l+z),/l+( z fi m ) + 



(l + z) 2 



-0. 



equal to 1.042 for our adopted redshift and cosmology. The 
concentration of the NFW halo, c, is then obtained as r^lrg. 

To provide better constraints on the mass distribution in 
the core of A478, we perform a finer binned "onion-peeling" 
analysis using just the Chandra data, which has much higher 
spatial resolution compared to XMM-Newton. The resulting 
gas density and temperature profiles are plotted in Figs. I10I& 
II II together with the best-fit 3rd order polynomial curve (in 
log-log space) to the data. These smoothed curves were used 
to evaluate the logarithmic gradients needed to determine the 
mass profile using equation^ Since we have fitted these poly- 
nomials in log-log space, the logarithmic gradient is simply 
the derivative of the function. 

The cumulative mass distribution was calculated at the ra- 
dius of the input T(r) and p(r) data, but using the values of the 
best-fit curve, rather than the measured data points. We have 
excluded the innermost bin since the gradient at this point is 
rather uncertain; we are more confident about the outermost 
bin since the turn-over implied by the best-fit curve is consis- 
tent with the temperature profile from XMM-Newton, which 
extends beyond this radius (see the diamonds in Fig.|6]i, and 
is also confirmed by the Chandra analysis of Vikhlin in et alJ 
(2004), which extends beyond the S3 chip. The mass profile 
is plotted in Fig.^J together with the best-fit mass function 
obtained with equation [5] The errors were determined from 
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FIG. 10. — High-resolution Chandra gas density profile. The solid line is 
the best-fit 3rd order polynomial relation in log-log space. 
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FIG. 1 1. — High-resolution Chandra gas temperature profile. The dashed 
and solid lines are the best-fit 2nd and 3rd order polynomial relations in log- 
log space, respectively. 



a series of 2000 Monte Carlo realizations of the input T(r) 
and p(r) data; in each case a best-fit mass model (equation[5j 
was fitted and the errors on the model parameters were ob- 
tained from the standard deviation of the set of points. Also 
plotted is the gas fraction profile (Fig. 1 1 31 . which actually 
decreases with radius out to ^200 kpc, beyond the confines 
of the central galaxy, somewhat contrary to expectati on (e.g. 
iDavid. Jones. & Formarl 19951 ISanderson et al.l20 03 ) . 

We find a best fit index parameter of n = 0.35 ± 0.22, with 
a corresponding scale radius and concentration parameter of 
r s = 317 ± 82 kpc and c = 6.9 ± 1.6, giving a v alue of R?m 
= 219 ± 1 25 kpc and R ? 5oa = 620 ±70 kp c. lAllen etalJ 
( 120031) and ISchmidt. Alle n. & Fabian (2004) also measured 
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FIG. 12. — Mass profile from the Chandra data. The solid line is the 
best fit using Equation|5] with an inner logarithmic slope of -0.35 ± 0.22; for 
comparison, the dotted line shows the best-fit NFW profile (inner logarithmic 
slope of -1). The \a errors were generated from a set of 2000 Monte Carlo 
simulations. 



the Chandra mass profile for A478, but fitted it with an 
NFW profile (i.e. n fixed at 1), to give R 2 oo = 2464+}^ 
kpc (converted to Hp=70) and 1?200 = 2367 kpc, respectively. 
Pointe couteau et alJ <2004l) also fitted an NFW profile to the 
poorer-resolution XMM-Newton data, to derive a value of /?20o 
= 2100 ± 100 kpc, in reasonable agreement with our findings. 
We estimate a total mass within ^?.n n of (1-3 ± .27) x 10 15 
M Q , compar ed to (1.84!°;^) M for lAllen et all and (1.1 x 
10 15 ) M for lPointecouteau et ail 

To investigate the robustness of our M(r) inner logarithmic 
slope measurement, we have also computed the mass profile 
for the case where a 2nd order polynomial is fitted to the tem- 
perature points (the dashed line in Fig. II It . In this case we 
recover a slightly steeper index of n = 0.49 ± 0.24. How- 
ever, we note that a 2nd order polynomial indicates a rising 
T(r) beyond the data points, whereas the 3rd order polyno- 
mial curve exhibits a turn-over that is clearly preferred by our 
XMM-Newton data, as well as the Chandra A CIS-I projected 
temperature profile of Yikhl inin et al.l ((2004). We return to 
discuss the mass distribution in ; 



6. DISCUSSION 
6.1. Hot spots in the ICM 

While there is generally a surprising amount of tempera- 
ture structure present in the ICM in A478, of particular note 
are four "hot spots" seen in the inner regions of the core (see 
Fig-0- Three of them are located to the North of the core; the 
other is to the South, and they all lie within ^30-60 kpc of the 
center. There are no known cluster galaxies located at the po- 
sition of these features, and they all lie just outside the edge 
of the central cluster galaxy, PGC014685: its d 2 5 diameter 
(the isophote at which the optical B band surface brightness 
reaches 25 mag/arcsec 2 ) is 57 kpc (Paturel et al. 1997). 

All four regions are clear outliers in both the temperature 
(Fig ■ IS} ar, d pressure (Fig.1141 profiles, but occupy a relatively 
small volume. Correspondingly they produce an essentially 
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FIG. 13. — Chandra gas fraction profile, which shows an unusual decrease 
with radius within the cool core, beyond the central galaxy. The \a errors 
were generated from a set of 2000 Monte Carlo simulations. 
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FIG. 14. — Pressure profile from the spectral mapping deprojection, with a 
smoothed curve shown as a solid line (see text for details). The hot spots are 
plotted as open circles. The X axis error bars have been omitted for clarity. 



negligible distortion in the X-ray surface brightness distribu- 
tion of the core of A478. The entropy of the hot spots also 
appears to be systematically higher than expected. Compared 
to the smoothed regression relation plotted in Fig.|9| there is 
an excess entropy of between 20-50 keV cm 2 present in each 
of the hot spots, corresponding to 30-50% of the expected 
value. The combination of excess pressure and entropy indi- 
cates that these are strong shocks, and therefore unlikely to be 
the result of sound waves or mildly supersonic motion. How- 
ever, the volume enclosing the hot spots is too small for them 
to produce a noticeable surface brightness contrast compared 
to the bright emission from the cool core. 

Some key properties of these hot spots are summarized 
in Table |4] including the estimated excess thermal energy 
present in each one, based on the excess gas pressure and es- 
timated volume. This excess was measured relative to an esti- 
mate of the underlying P(r), obtained from a locally weighted 
regression of the data in log-log space (the solid curve in 
Fig. I14i. using the method described in %5.3\ The excess en- 
ergy in each hot spot was calculated by multiplying its resid- 
ual excess pressure above this curve by the estimated volume 
for the region. Although the errors on the individual estimates 
are quite large, the total excess thermal energy contained in all 
4 hot spots combined is (3 ± 1) x 10 59 ergs. We note that this 
result depends only weakly on the volume, V, assumed for the 
hot spots: the energy, E oc V p, and p oc y/K/V (Equation^, 
where K is the MEKAL model normalization, which implies 
E oc \fKV. 

The global morphology of A478 exhibits a regular struc- 
ture, and the XMM-Newton hardness temperature map (top 
right panel of Fig.Q shows no sign of significant disturbance 
in the outskirts of the ICM. This, combined with the presence 
of a well established cool core, suggests that A478 cannot 
have experienced a recent merger or significant disturbance, 
which might account for the existence of the hot spots. It can 
be seen from Table|4]that the distances from the cluster center 
and energy excesses are similar for all the hot spots. Further- 



more, they appear to be located (see top left panel of Fig. [7} 
in rough alignment with the (approximately) conical depres- 
sions discovered bv lS03l corresponding to radio lobes emitted 
from the active nucleus in the central galaxy. These charac- 
teristics point to the influence of an event originating in the 
cluster center. Correspondingly, it is reasonable to associate 
them with AGN activity. 

On the basis that the hot spots share a common origin, we 
have extracted a combined spectrum from all 4 regions and 
fitted it with a single absorbed MEKAL model. In this way 
we can gain a better understanding of the degree to which 
the gas has been heated. We obtain a best fit temperature of 
12.1*2 6 ^ e ^ ; an abundance of 0.26*Q26 solar ar, d aborption of 
(2 .4*q 2 x 10 21 ) cm" 2 . However, by fixing the temperature to 
the ambient value of 5.5 keV (e.g. see the best-fit curve at ~ 
50 kpc in Fig. II 1> and refitting with the other parameters left 
free to vary, we obtain a change in fit statistic of A\ 2 = 16.7, 
corresponding to 4.1cr significance. 

We stress that the interpretation of this result as evidence 
for the existence of the hot spots hinges on the validity of 
associating these 4 separate regions with each other. It is cer- 
tainly possible to achieve a similar level of significance by 
randomly grouping hot regions arising from purely statisti- 
cal fluctuations. However, we believe that the location of the 
hot spots in the cluster core points to an association between 
them that indicates they are not merely random peaks in the 
temperature distribution. Specifically, they lie within a narrow 
annulus centered on the cluster peak, and they are positioned 
within the range of angles that enclose the roughly conical bi- 
polar cavities excavated by radio lobes from the central AGN. 

6.2. AGN Heating? 

If the hot spots are indeed the result of AGN activity, we can 
address the issue of how they may be connected to the cavi- 
ties seen in the ICM. The fact that the hot spots are located at 
least 3-4 times further away from the cluster center implies 
a different formation mechanism to that which excavated the 
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TABLE 4 

Data for the hot spots in the central ICM of A478 



Region 


Distance" (kpc) 


Cts 6 


kT (keV) 


X 2 /dof 


L x c (10 42 erg/s) 


Vol rf (10 68 cm 3 ) 


AE e (10 58 erg) 


'coo/ (yr) 


A 


48-59 


1054 


|2 1+2.6 


57.65/43 


9.3 


2.3 


11 ± 6 


6xl0 8 


B 


34-42 


727 




28.93/29 


6.7 


0.9 


4.6±4 


4xl0 8 


C 


38-51 


1108 




45.08/45 


10 


2.6 


8.7 ±6.5 


6xl0 8 


D 


45-55 


639 


10.3:*! 


25.99/25 


5.4 


1.2 


5.0±4.6 


5xl0 8 



NOTE. — "Nearest/furthest distance to hotspot, measured from the cluster center; *net source counts in the 0.5-7.0 keV band; 
c bolometric X-ray luminosity; ^estimated volume of region; c estimated excess energy in hot spot; -'gas cooling time (E/Lx, using 
the total rather than the excess hot spot energy). All errors are la. 



cavities. In addition, S03 estimate that the minimum energy 
needed to create these cavities is ~3x 10 58 ergs ( Birz anet al.l 
( 2004) estimate ~ 1 .2 x 10 58 ergs), which is a tenth of the ther- 
mal energy contained in the hot spots. IS03I further estimate 
a cavity age of ^3x 10 7 yr, from the time taken for the bub- 
bles which formed t hem to rise buoyantly to their observed 
location. Similarly, Birz an et"aI1 estimate a cavity age of be- 
tween 1-3 x 10 7 yr, by considering 3 different timescales. For 
comparison, the time taken to reach a distance of 50 kpc at the 
local sound speed (920 km s" 1 ) is ~5 x 10 7 yr, for a mean am- 
bient gas temp erature of 5 keV (c s =1480(r g /10 8 K) 1 ' 2 km s" 1 
(Sarazin 1988)). There is no evidence of a bow shock or X- 
ray emission excess anywhere near these features, which ar- 
gues against very recent supersonic motion. It is more likely 
that the hot spots may have been generated in situ, perhaps in 
a similar manner to radio hot spots located at the shock ter- 
mination of radio jets, which have been observed at similar 
distan ces from the cluster cen ter (e.g. see the recent sample of 
IHardcastle & Sakellioul2004l) . 

An estimate of the minimum time scale for heat conduc- 
tion to erase the hot spots can be obtained by considering the 
saturated heat flux, given by 

/2kT \ 1/2 

<?sat = 0.4 n e kT e (6) 

V TTJTCe / 

dCowie & McKedl 19771) . where m e is the electron mass. For 
hotspot A the gas density is 2 x 10 -4 cm 3 and its temperature 
is 12 keV. Assuming a spherical geometry for this region, with 
a volume of 2.3 x 10 68 cm 3 (Table0}, yields a cooling rate of 
10 44 erg s" 1 . Thus, it would take ~ 3 x 10 7 yr to lose its ex- 
cess energy, of 10 59 erg, via conduction. This is about a tenth 
of the likely lifetime of the hot spots, which implies that heat 
conduction must be suppressed by a factor of ~10 compared 
to the Spitzer rate. Such a decrease in conduction efficiency 
could easily be achieved by magnetic fields, which have been 
shown to lead to suppressi on factors of ~5 in a w eakly colli- 
sional magnetized plasma (Naravan & Medvedev 2001). 

The positions of the Northern set of 3 hot spots (Fig. [7} 
may offer a clue as to the mechanism responsible for their 
creation. Their spacing is suggestive of a wide jet opening 
angle, of roughly 100 deg rees. This is consistent w ith the 
slow jet heating scenario of Soker & Pizzolato (2004), which 
relies on such a wide angle to produce a poorly collimated, 
massive and relatively slow outflow. Moreover, the emissiv- 
ity of the small bubbles in the ICM predicted from this model 
is too low for them to be detected in X-ray images with ex- 
isting telescopes, c onsistent with the hot spots in A478. The 
ISoker & Pizzolato m odel is similar to the effervescent heating 
model of lBTgehnanl (120011) . which also postulates the forma- 



tion of many small bubbles, which heat the ICM by doing 
PdV work as they expand and rise. A particularly attractive 
aspect of this scenario is that the heating is smoothly regu- 
lated, since gas cooling acts to steepen the ICM pressure pro- 
file, thus yielding more energy from the bubbles' expansion 
so as to balance the cooling. However, this mechanism for 
heating the ICM would need to be slow enough to avoid dis- 
rupting the nearly power-law scaling of entropy with radius 
seen in Fig. |9] 

While the interaction of AGN with the intracluster medium 
has b e en well studied in recent simulations ( e.g.lRizza et al. l 
2000; Reynolds, Heinz, & Begelman 2001; Churazov et al. 
200 11 Briiggen & Kaiserl 120021: iBasson & Alexander] 120031 
Ruszkows ki et alJ200 4a). the formation of shock heated blobs 
of the type seen in A478 does not appear t o be a generic 
feature. However, Clarke, Harris, & Carillt i ll 9971) present 
3D magnetohydrodynamical simulations of a supersonic jet, 
which predict the formation of a pair of X-ray excesses, cor- 
respondingly to shocked ambient gas, at either side of the jet 
orifice, which itself produces an X-ray cavity of the type seen 
in A478. The iClarke et al.l simulations proceed for ^10 7 yr, 
which does not allow enough time to track the full evolution 
of these excesses. However, it seems probable that the shock 
heated gas in these blobs would expand outwards and heat 
the ICM, which could be an important means of transferring 
AGN mechanical luminosity into gas thermal en ergy. 

Sim ilarly, the more recent simulations of lOmma et alJ 
(2004) appear to produce three distinct density enhancements 
surrounding the head of an AGN outflow, appearing after 
~ 10 8 yr. However, the jet responsible for these features is 
only weakly relativistic, so may not be capable of heating 
these blobs at the level observed in A478. Observationally, 
perhaps the closest analogue to the hot spot s reported here 
is the heated bubble observed in MKW 3s by Ma zzotta et alJ 
(2002). This feature lies at roughly 90 kpc from the cluster 
core and is approximatel y 25 kpc in radi us in the plane of the 
sky. More recently, Maz zotta et alJ ( 120041) present new, lower 
frequency radio observations for this cluster, which show that 
the radio emission is fully enclosed by the shocked gas. A 
similar radio/X-ray interaction is observed in Hydra A, where 
expanding radio lobe s, observed at 330 MHz, appear to be 
driving a shock front (Nulsen et al. 2004J). It is possible that 
equivalent lower frequency observations of A478 may reveal 
larger radi o lobes than tho se seen in the 1.4 GHz VLA map 
(shown in lSun et alJ2 003). which might show some evidence 
of an interaction with the hotspot regions. 

However, unlike A478, the MKW 3S bubble is clearly vis- 
ible as a depression in the surface brightness and the Hydra 
A shock front is revealed as a discontinuity in the X-ray im- 
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age. Similarly, there are spiral arm-like features observed 
by Chandra in the elliptical galaxy NGC 4636, produced by 
weak shocks driven by off center outbursts JJones et all 2002). 
However, these are ostensibly density enhancements, which 
produce a clear signature in the surface brightness - as already 
noted, the hot spots seen in A478 are effectively invisible in 
X-ray images. 

There is further evidence of possible AGN activity disturb- 
ing the gas in the form of a tentative "ripple" in the surface 
brightness profile of the core (Fig. ll5l l. We have generated an 
unsharp mask image, using the method of Rusz kowski et alJ 
(2004b), who have simulated viscous dissipation in the ICM 
caused by AGN activity. A 0.5-7.0 keV raw image was 
smoothed with a Gaussian kernel of FWHM 6 kpc (a=3.1 
pixels) and subtracted from the unsmoothed image; the re- 
sulting image has been smoothed with a Gaussian kernel of 
cr=1.5 pixels to highlight the features present. The depres- 
sions caused by the radio lobes are clearly seen, producing an 
"H" sh ape which resembles that seen in the Ruszkowski et al. 
(2004b) simulations, and there is also a suggestion of a "rip- 
ple" to the North of the center, lying between 16 and 26 kpc 
of the AGN. This ripple has a low contrast, but is clearly 
longer and more coherent that the typical noise features in 
the image. I t is reminiscent of the ripples seen in the ICM of 
the Perseus ( Fabian et al. 2003) and Virgo (at 14 and 17 kpc) 
llForman et al J 12004) clusters, which are attributed to shock 
fronts emanating from the AGN in the central cluster galaxy. 
If this ripple is genuine, the implied travel time for the feature 
is roughly 2x 10 7 yr (given the assumed local sound speed of 
920 km s" 1 used above), which is comparable to the age of the 
cavities. We note that the hot spots in A478 are not visible in 
the unsharp mask image and that they lie at least twice as far 
out as the putative ripple. 

6.3. Entropy Profile 

The radial variation in gas entropy in A478 (Fig. [9} is un- 
usual, in that it appears to follow a broken power law rela- 
tion, with an outer logarithmic slope of ~1.1 and a slightly 
shallower inner slope of ^0.95 within ^0.1 R500, but with no 
indication of a flattening in the core. Non cool-core clusters 
generally exhibit a flatte ning of the entropy profile towards 
the ce nter (e.g. [Rasmussen & Ponman 2004; iPratt & Arnaudl 
2005), while the effect of cooling is to lower the entropy in 
exactly this region. It therefore appears to be a conspiracy 
that causes these two effects to counterbalance each other 
so as to preserve the scaling evident in the outskirts, where 
cooling is insignificant. A recent XMM-Newton study of 13 
nearby cool core clusters has found that their entropy profiles 
follow a power-law with best-fit logarithmi c slope of 0.95 , 
out to roughly half the virial radius (Piffa retti et alJl2005l) . 
However, other cool core clusters show evidence of flatten- 
ing in their core entropy profile (e.g. I Pratt & Arnaudll2003t 
lO'Sullivan et all2003tfwise et all2004l) . Moreover, the slope 
of ~ 1 . 1 is p redicted to arise from shock heating due t o accre- 
tion of gas (Toz zi et alJ l2000: Tozz i & Normanl l200l . which 
would not be expected to persist in cooling-dominated cluster 
cores. 

The shape of A478's entropy profile resembles that of the 
fossil group NGC 6482 (kT ~ 0.7 keV), which was the sub - 
ject of a recent Chandra analysis ( iKho sroshahi et alJl2004l) . 
That observation is restricted to the inner 10% of Z?200 of the 
halo, but is is clear that the entropy data are consistent with 
a power law form, with a logarithmic slope slightly shallower 
than S oc r L1 . Since NGC 6482 is a fossil group, it is expected 




FIG. 15. — An unsharp masked 0.5-7.0 keV Chandra image of the core of 
A478. Each pixel is 0.492" (0.8 kpc) across and the image has been smoothed 
with a Gaussian of a = 1.5 pixels. The edges around the radio lobe cavities 
are clearly visible as an "H" shape, and there is evidence of a possible curved 
"ripple" to the North of the cavities. The positions of the other features are 
labeled exactly as positioned in Fig. 171 



to be a very old system. However, it shows no evidence for 
strong cooling and, furthermore, has a gas temperature pro- 
file which decreases monotonically with radius, thus ruling 
out thermal conduction as a means of heating the core. Since 
NGC 6482 shows some evidence of having a mildly active 
nucleus, it is possible that this system too has been subjected 
to a phase of AGN heating which has only recently abated, 
where significant gas cooling is presently being reestablished. 

Given the apparently significant impact of the AGN on the 
gas in the core of A478, it is worth noting that excessively 
large heat input would trigger convective motions that would 
stir up the gas, flattening the entrop y profile and erasing abun- 
dance gradients in the core dChurazov et alJl2001t [Briiggen 
2002; IChurazov et al. As there is some enhancement 

in the metallicity towards the center the ICM (Fig. |2}, it is 
unlikely that a significant disturbance could have occurred re- 
cently. This may be because the AGN heating is widely dis- 
tributed (e.g. effervescent) and thus not very vigorous. Alter- 
natively, it could indicate that the most unstable (i.e. lowest 
entropy) gas is also the most metal rich, and that as it cools 
and moves inwards it produces a central abundance peak. 

The cooling time of gas in the center of A478 is only 
~9 x 1 7 yr — roughly 6 times shorter than the cooling time of 
the hot spots (see Table0). S03 point out that the observed ra- 
dio lobes are much smaller than the X-ray cavities, which sug- 
gests a fading of the radio source on a timescale of ^10 7-8 yr. 
Furthermore, the ratio of mechanical power needed to exca- 
vate th e cavities to the cu rrent radio power of the AGN is over 
1000 (iBirzan et al.l2004l) - by far the largest ratio amongst the 
10 clusters (from their sample of 18 objects) for which they 
have such data. The inclusion of the hot spot contribution to 
the energy budget raises this value even higher. Therefore, it 
is possible that we are witnessing the rapid reestablishment of 
a cluster cooling flow, following a period of AGN heating of 
the gas, before the cavities have been completely refilled. On 
the other hand, an AGN can vary in radio power w hile its jet 
power remains constant (Eilek 2004; Birzan et al. 2004), and 
thus some mechanical heating of the ICM may be ongoing, 
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despite the fading radio luminosity, as possibly indicated by 
the putative "ripple" seen in Fig. 1151 

If similar hot spot features are present in other clusters with 
radio lobe cavities, the timescale for which these features per- 
sist as over-dense blobs is short enough to make them unlikely 
to be detected in X-ray images. Furthermore, for much of this 
time they would reside in close proximity to the cluster cen- 
ter, where emission from the cluster core would be substan- 
tial, greatly reducing their surface brightness contrast. The 
alternative method of detecting them - locating the hot spots 
from their temperature signature - is also challenging, since it 
requires high resolution hardness mapping and the best pos- 
sible data quality (the A478 Chandra observation analysed in 
this work comprises ^400,000 net source counts from the S3 
chip). 

6.4. Mass Profile 

The soft core in the mass profile (Fig. I12> of A478 is 
unusual, and apparently at odds with numerical simulations 
of clusters, which i ndicate inner logarithmic slopes of -1 
llNavarro et al] [l995) or even as steep as -1.5 ( Moo re et alJ 



1998). However, such a flat slope is permitted by purely ana- 
lytical considerations of t he dark matter , base d on fundamen- 
tal statistical mechanics: lHansen et alJ ll2005l) have recently 
demonstrated that this approach permits < n < 10/3, for 
the index parameter in equation |5] Furthermore, a number 
of clusters have b een found with similarly flatter cusps in 
their mass profiles. Tvson. Kochanski. & dell' Antonio ( 199 
measu red a slope of 0.57 ± 0.02, and ISand. Treu. & Elli; 
(2002) observed a value of 0.35 for two different clusters, 
based on strong gravitational lensing analyses. More recently, 
ISand et a"fl ll2004l) have studied lensing arcs in 6 clusters and 
report a mean slope of 0.52 ±0.05. It is interesting to note that 
their sample comprises only clusters with a dominant bright- 
est cluster galaxy (BCG), as is the case for A478. Also, the 
cluster Abell 1795 has been shown to have an inner mass pro- 
file logarithmic slope of 0.59*° }t (90% con fidence), inferred 
from a Chandra analysis (Ett ori et al.l20 02). 

It must be remembered that, unlike lensing measurements, 
our X-ray analysis hinges on the assumption of hydrostatic 
equilibrium, which could easily break down in the core, 
given the impact of the central AGN. For example, the pres- 
ence of significant non-thermal pressure support could mimic 
a less concentrated mass distribution. The simulations of 
iFaltenbacher et all J2005) indicate that random bulk motions 
of gas account for 10% of the total pressure support in clus- 
ters. Also, a recent XMM-Newton analysis of the spectrum 
of pressure fluctuations in the Coma Cluster found a lower 
limit of ^10% for the contribution to the ICM pressure sup- 
port arising from turbulence. It is possible that bulk or tur- 
bulent motion of gas associated with AGN activity could be 
significant in the core of A478. Presumably such a mecha- 
nism would give rise to a radially diminishing pressure con- 
tribution, which could act to flatten the inferred mass distribu- 
tion in the manner observed. In such a situation, the inferred 
gas fraction would be increasingly overestimated at smaller 
radii, which may explain the radially decreasing trend ob- 



served within ^200 kpc in Fig. 1131 

7. CONCLUSIONS 

We have studied the detailed thermodynamic properties 
of the intracluster medium in the relaxed, cool-core cluster 
Abell 478, using Chandra and XMM-Newton X-ray data. Our 
main findings can be summarized as follows: 

1 . We find that our Chandra X-ray temperature measure- 
ments are systematically h otter than those from XMM- 
Newton, as also reported by Pointecouteau et al. ( 2004). 
By selecting approximately isothermal regions, we find 
slightly better agreement between the two. However, 
by simulating multiphase spectra and fitting them with 
a single temperature model, we find no evidence for 
significant disagreement between Chandra and XMM- 
Newton. We therefore conclude that the observed dis- 
crepancy cannot be fully attributed to non-isothermality 
in fitting single-temperature models. 

2. The entropy profile appears to agree well with the 
empiri cal modified entropy scaling of iPonman et alJ 
(2003). Moreover, the power law trend continues to 
the innermost radius measured (<10 kpc), with only a 
slightly shallower slope. There is no evidence of any 
core in the entropy profile. 

3. Under the assumption of hydrostatic equilibrium, we 
infer a mass profile that exhibits a soft core, character- 
ized by a logarithmic slope of-0.35±0.22. This is sig- 
nificantly flatter than an NFW profile (slope = -1), but 
is consistent w ith the recent gravitational lensing results 
of Sand et al. (2004) for clusters containing a dominant 
central galaxy. 

4. We have discovered four hot spots in the ICM located 
well within the cool core, where the gas is roughly 
twice as hot as its surroundings. The combined excess 
energy associated with these regions is (3± 1) x 10 59 
erg, which is ^10 times the energy needed to exca- 
vate the cavities produced by radio lobes from the cen- 
tral AGN. The properties of these hot spots suggest 
they may the result of strong shock heating from a 
jet/outflow originating in the AGN. 
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